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The Gauss-Bonnet (GB) curvature invariant coupled to a scalar field (j) can lead to an exit from a 
scaling matter-dominated epoch to a late-time accelerated expansion, which is attractive to alleviate 
the coincident problem of dark energy. We derive the condition for the existence of cosmological 
scaling solutions in the presence of the GB coupling for a general scalar-field Lagrangian density 
p{4>,X), where X = — (1/2) {\7(j>)^ is a kinematic term of the scalar field. The GB coupling and the 
, Lagrangian density are restricted to be in the form f{4>) oc e^** and p — Xg{Xe^'^), respectively, 

■ where A is a constant and g is an arbitrary function. We also derive fixed points for such a scaling 

' Lagrangian with a GB coupling ficj)) oc e^** and clarify the conditions under which the scaling 

CN) ^ matter era is followed by a de-Sitter solution which can appear in the presence of the GB coupling. 

^ . Among scaling models proposed in the current literature, we find that the models which allow such 

' a cosmological evolution are an ordinary scalar field with an exponential potential and a tachyon 
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field with an inverse square potential, although the latter requires a coupling between dark energy 
and dark matter. 

PACS numbers; 98.70.Vc 
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; I. INTRODUCTION 



A great deal of efforts in modern cosmology, over the course of the past decade, have been made in trying to 
understand the nature and the origin of dark energy vifhich is responsible for late-time acceleration of the universe 
! PjHjUI- Although cosmological constant is the simplest candidate of dark energy consistent with current observations, 
it suffers from a severe fine-tuning problem from a field theoretic point of view [j] . It is then tempting to find alternative 
'"pi models of dark energy which have dynamical nature unlike cosmological constant. An accelerated expansion can be 
' realized by using a scalar field whose origin may be found in superstring or supergravity theories. Thus modern particle 
physics has been an attractive supplier for the construction of scalar-field dark energy models-such as quintessence 
(— I I Q , K-essence @ , tachyon Q , dilatonic models Q , and ghost condensate d, • 

■ Although the dynamically evolving field models have an edge over the cosmological constant scenario, they too, in 
general, involve fine-tuning of initial values of the field and the parameters of the model. Secondly, nucleosynthesis 
places stringent restriction on any additional degree of freedom (over and above the particle physics standard model 
5-H ' degrees of freedom) which translates into a constraint on the ratio of the field energy density to the energy density 
of the background fiuid 3]. These problems can be addressed by employing scalar field models exhibiting scaling 
solutions [ill . IT3 . [Tp. [tS , [isl . [l6j . For instance, a standard scalar field with an exponential potential in General 
Relativity (GR) [III, can mimic the background (radiation/matter) remaining sub-dominant so as to respect the 
nucleosynthesis constraint [l3|- The scaling solutions as dynamical attractors can successfully address the second 
problem and can also considerably alleviate the fine-tuning problem of initial conditions. 

For a viable cosmological evolution, the scalar field should remain unimportant for most of the cosmic evolution 
and should emerge from the hiding epoch only at late times to give rise to the accelerated expansion. Since scaling 
solutions are non-accelerating, the model should be supplemented by an additional mechanism allowing the field to 
exit from the scaling regime at late times. There are several ways of achieving this goal: 

• (i) At late times the potential becomes shallow to lead to the accelerated expansion. For example this is realized 
by the double exponential potential given by V{(p) = Vo{er^'"l' + e"'""^) with > 20 and < 2 [l3| (here 
= SttG with G being gravitational constant). 
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• (ii) A local minimum of the field potential appears at late times. The examples of this type of models are 
V{(f>) = Vb[cosh(AK<?!)) - 1]" [ISi] and V{<j)) = Foe"^'"^[A + - B)^] In the former case the acceleration 
occurs for n < 1/2 on average, whereas in the latter case the system always approaches a de-Sitter universe 
after a damped oscillation of the field. 

• (iii) Assisted quintessence 120j. This is realized by considering many fields with the sum of the exponential 
potential V = V A^e"'*'''*''^', since the presence of many scalar fields leads to the smaller effective coupling, i.e.. 

Although these are interesting attempts to provide a successful exit from the scaling regime, it is fair to say that 
the aforementioned mechanisms are still phenomcnological. It is of great importance to find cosmological models from 
fundamental theories of particle physics like string theory which may incorporate the findings of the phenomcnological 
studies. String theory necessarily includes scalar fields (dilaton & modulus fields); the low energy string effective 
action also gives rise to a modification of standard GR in terms of Riemann invariants coupled to scalar field (s) 
[2^ . The modifications are important to address the problems of past [2^ and future [l^] singularities. Out of these 
corrections, the Gauss-Bonnet (GB) term is of a special relevance. The GB term is a topological invariant quantity 
which contributes to the dynamics in 4 dimensions provided it is coupled to a dynamically evolving scalar field. It 
is the unique invariant for which the highest (second) derivative occurs linearly in the equations of motion thereby 
ensuring the uniqueness of solutions. Nojiri et al. [2a | studied scalar-Gauss-Bonnet cosmology to demonstrate the 
existence of a particular dark energy solution. Such a solution exists in the presence of an exponential potential 
V{(j)) oc e^'^'"^ and the GB coupling /(^) oc e^^'f' . More recently Koivisto and Mota [1^ showed that the GB coupling 
allows a possibility to lead to the exit from the scaling matter era to the dark energy dominated universe for the same 
model. See Refs. ^27i, i28, 2^, jiOj, HI, HI, [H, |35., JBg,] for other aspects of the GB dark energy scenario. 

In this paper we shall clarify the form of the GB coupling /(0) for the existence of scaling solutions by starting from 
a very general scalar-field Lagrangian density p{4>, X) where X is a kinematic term of the field. This includes a wide 
variety of scalar-field models such as quintessence, K-essence, tachyon, and ghost condensate. The demand for scaling 
solutions not only determines the field Lagrangian to be X) — Xg{Xe^'^) [g is an arbitrary function) but also 
fixes the GB coupling in the form /((/)) (x C^*^ (here we use the unit = 1). Thus our results are useful to construct 
viable dark energy models which possess scaling solutions in the matter-dominated epoch. We then investigate the 
fixed points for the scaling Lagrangian p((/), X) = Xg{Xe^'^) with a GB coupling /(</>) cx e'"^. We show the existence 
of a pure de-Sitter fixed point which exists only for specific models and study a possibility to obtain a scaling matter 
era that finally approaches the de-Sitter solution. A standard field with an exponential potential is a specific model 
to realize such a cosmological evolution provided that /i > A. A tachyon field also possesses the de-Sitter fixed point, 
but it does not have a scaling matter epoch. However the scaling matter era followed by the de-Sitter solution is 
possible if dark energy couples to dark matter. We shall discuss cosmological dynamics in such cases in details. 



II. THE CONDITION FOR THE EXISTENCE OF SCALING SOLUTIONS IN GAUSS-BONNET 

COSMOLOGY 



The model we study is given by the following action 



S = d Xy/-gM 



-LR + p^X,cl>)-f{^)Rl^ 



(1) 



where gM is a metric determinant, i? is a Ricci scalar, ^(0, X) is a Lagrangian density in terms of a scalar field and 
a kinematic energy X = -{l/2)g'"'df,(j) d^cj). The Gauss-Bonnet (GB) term, i?^B = i?^ _ iR^^Rt^'' + R^^f.^R'^l^f" , 
couples to the field (f> through the coupling f{(j>). We allow for an arbitrary coupling between the matter fields ipi and 
the scalar field (j). We assume that the field (p is coupled to a barotropic perfect fiuid (energy density pm and pressure 
density Pm) with a coupling given by (37i. i38>. i39i] 



1 



6S,, 



Pm\/ — 9M 



(2) 



Note that the energy density of the scalar field is given by p = 2X{dp/dX) — p with an equation of state = p/p. 
In what follows we shall use the unit = SirG = 1. 

In the flat Friedmann-Robertson- Walker background with a scale factor a, the equations for p and p^ are 



p + 3H{1 + w^)p = -Q/9,„0 - 24/'(0)0i/2(ij2 ^ ^) 

+ 3i?(l + Wm)Prn = +QPrn'i' , 



(3) 
(4) 
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where H = a/a is a Hubble parameter. A dot and a prime represent derivatives in terms of cosmic time t and the 
field (j), respectively. The Hubble parameter satisfies the constraint equation: 



(5) 



Then the energy density of the GB term is given by pcB = '^^f {4>)'i>H^ ■ We define the fractional densities of the field 
0, the GB term, and the matter fluid, respectively, as 



GB 



4i — "GB ■ 



(6) 



Scaling solutions arc characterised by constant values of O^, Ogb and f2m, or equivalently a constant ratio between 
each energy density. Then this gives the condition: 



d In p d In p„i d In pcB 



-3(f + Woft) , 



dA^ diV dN 
where WcS is an effective equation of state related to the Hubble rate via the relation 

H 3, 



(7) 



(8) 



We also note that and Wm are assumed to be constants in the scaling regime. From Eqs. ([3]) and (jlj together with 
the scaling condition ([7]), we find 



^d^ 



^0 ~1~ 

Substituting this relation for Eq. (|3|) and using the relations ([7]) and ([S]), we obtain 



Wcff = 



1 + ri^ + Vtm 



Hence Eq. ([9]) can be rewritten as 

f 



■ diV + 17„ 



3(wm - W^)Q.4, + Q. 



GB- 



Vt^{\ + 3w0) + Sl,„(l + 3w„i) 



From these equations and the definition of X, we find 



which gives 



2X = 



A\nX 



Q 



1 + Jl^ + Vtra 

p{X,<f) 



= C {— const) . 



dlnQ 



dN ' ' dN 

From Eq. ([7]) together with p oc p, the Lagrangian density p — p{(t>, X) satisfies the differential equation 

d\np dX d In p d( 



dX dN d<j> dN 
Substituting Eqs. (ITT|) and for Eq. IH]), we arrive at 



-3(1 + ?z;eff) 



1 



2 dQ(0) 
\Q d(t) 



dlnp 1 dlnp 
dlnX ^ A dch 



where 



2(1 - »gb)[3(1 + Wm^^ + w^n^) - 2»gb] 

3{Wm - W^)fl^{2 - JIgb) + f^GB[l - f^GB + 3(w„ri,„ + IW^ri^)] 



(9) 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 



Q(0). (16) 
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Equation ([15]) is the same form of equation which was obtained in the absence of the GB coupling ^Td\ . Following the 
procedure in Ref. fld\ , Eq. restricts the form of the Lagrangian density to be 



(17) 



where g is an arbitrary function. If Q is a constant, this reduces to a simple form: 

p = Xg{Y), r = Xe^^ (18) 

where p is redefined so that a constant factor is removed from Eq. p7l) . This form of the scaling Lagrangian was first 
derived in Ref. [l3| in the absence of the GB coupling (see also Ref. [ij] for the extension of the analysis to the case 
(X p"). Even if the GB term is present, we have shown that the Lagrangian takes the same form with a modified 
value of A given in Eq. (I16p . For a later convenience we present several models which belong to the scaling Lagrangian 

m- 

• (i) An ordinary field with an exponential potential, i.e., p = X — ce~^'^ [l3|. This corresponds to the choice 

g{Y) = l-c/Y. (19) 

• (ii) A dilatonic ghost condensate model, i.e., p = —X + ce^'^X'^ [l3|. This corresponds to the choice 

giY)^-l + cY. (20) 



(in) A tachyon field, i.e., p = —V{^p)^yl — tp'^ with V{ip) oup [40j EH l43| . This corresponds to the choice 



g{Y) = -{c/Y)VT^ . (21) 
Note that the field if is related to via the relation tp = (2/A)e'^'^/2. 

From Eq. (|11|) we find 4)/H — C/Q. Since the GB energy density, pcB — '^4:f'{(j))(pH^, satisfies the relation 
din pqb / <iN ~ — 3(1 + Wcs) for scaling solutions, one gets 

/''(0) on _L ^ 

-—-^3{l + u..^). (22) 

This yields f" {(f) / f {(p) = 3(1 + Wcff)i5/C = A, which gives 

/'(0) = ae^* , (23) 

where a is an integration constant. This means that the GB coupling is restricted to be in the form /(0) oc e^*^ for 
the existence of scaling solutions. 

In the case of a non-relativistic dark matter {wm ~ 0), the quantity C is given by 

^ _ 1 - - ^^r« - &w^Vt^ ^2^-) 

1 + fi^ + Vim 

whereas the effective equation of state is 

^ 3(l + 0, + 0™) ■ 

Hence one has a simple relation = — C/3. From Eq. (fTB]) we get C = 'SQ/{Q + A). Then from Eq. pT|) we find 
the following relation along the scaling solution: 

(26) 



H Q + X 

This is the same relation as in the case where the GB term is absent [l^, [13] • 
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III. AUTONOMOUS EQUATIONS AND FIXED POINTS 



In this section we shall obtain fixed points for the scaling Lagrangian (|T8)) . We take into account the contribution 
of both non-relativistic dark matter {w^ — 0) and radiation {wr — 1/3), in which case the Friedmann equation is 



(27) 



Here the energy density of the scalar field is given by p = X{g + gi), where g„ = with g'-'^^Y) being the 

n-th derivative in terms of Y. We assume that the field 4> is coupled to dark matter with a coupling Q and that it 
is uncoupled to radiation, as it is the case in the context of scalar-tensor theories [s^. Hence the equations for the 
energy densities of dark matter and radiation are 



Pm + iHpm = +Qpm4> , 
Pr + AHpr = . 

The scalar-field equation ([3|) can be written as 

'<}> + ?>AH{g + gi)4> + \X[\- A{g + 2.gi)] -f^ 2Af' {(t>)AH^ {H^ + H) + AQp„, = , 
where A is defined by 

A={g + 5gi + 252)"' • 
We shall define the following dimensionless quantities: 



Xi 



Then the energy fractions are 

^4, = xlig + 2gi) , r^GB = 8\/6a;i2;3 , flr=xl, il^ 
The variable Y defined in Eq. can be expressed as 

Y = X1/X2 ■ 



Xi 



1 — 17^ — JIgb 



(28) 
(29) 



(30) 
(31) 

(32) 

(33) 
(34) 



The GB coupling /'((/>) takes the form for the existence of scaling solutions. We shall derive autonomous equations 
for a more general case: 



f{4>) = ae^* . 
Since p — X for scaling solutions, one has the relation 

X2X3 = a/3 . 

In this case the variable X3 is not needed to study a dynamical system. 

Taking the time-derivative of Eq. ^7} with the use of Eqs. (|28l) and we obtain 



2H = -2X(g + gi) p,„ - -pr 



(35) 



(36) 



(37) 



Eliminating the (j) term from Eqs. ([30)l and ([37|) . one gets the fohowing relation 

H 1 

(1 - 8\/6a;ia;3 -|-96Aa;^)-^ = --{3 + 3gxl + xl) 



+4X3 



2V6X1 -I- 3(2^ - X)xl + 3a|-V6(5 + gi 



)xi + (Q + X)xl{g + 2gi) 



-8x3 + Q{xl -\) + BV&Qxi 



(38) 
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The autonomous equations are 

-3A{g + gi)xi 



dN 



dX2 

dN 

dxs 
dN 

dx4 
dN 



V6 
2 



A{Q + X)ig + 2gi)xl - Xxf + QA{xl - 1 + SVexiXg) 



'4:V6Ax3 - {xi + aVgAxs)— , 




(39) 
(40) 

(41) 

(42) 



where H/H"^ is given in Eq. ([38|) . 

In what follows we will obtain fixed points in the absence of radiation (xi = 0). Let us consider the cases: (i) A = /i 
and (ii) X ^ fi separately. We shall study the case A > without loss of generality. 



A. Case of A = /i 



As we showed in the previous section, there exist scaling solutions for A = /i. Taking note of the relation (|36p . 
neither X2 nor x^ is identical to zero if a 7^ 0. When a is negligibly small, X2 and X3 can be very small as well. In 
this case we may regard 0:2 ~ or 0:3 « as approximate fixed points in Eqs. ()40p and ()41|) . Here we do not consider 
such fixed points, but in the next subsection we will discuss those cases when A does not equal to /z. 

Then from Eqs. (|40)) and (|4T|) . the critical points corresponding to a;2 7^ and 0:3 7^ satisfy 



H 



V6 



Xxi 



Then substituting Eq. ^ for Eq. we get 

(g + X){g + 2gi)xl - Q - V&{g + gi)xi + 4x3 V6(2Q + X)xi - 2 
From Eqs. ^ and (gS]), we find 

?)gx\ — V^Xxi + 3 = %xiX3,{2V& — 3Axi) . 
We note that the following useful relations hold for the scaling Lagrangian (US]): 



= 0. 



g 



w^Q^ ^ gx^ , g + gi 



(1 + w^)Q^ 



Substituting these expressions for Eqs. and ([35]) and eliminating the w^^^ term, we obtain 

2(Q + A)a;i - VE] {1 - - Qgb) 



0. 



(43) 



(44) 



(45) 



(46) 



(47) 



This gives the following fixed points: 
• (a) Scaling solution: x- — — — 



2{Q+\) ' 

• (b) Scalar- field & GB dominated point: fl^ + Hqb = 1 • 
Since both points satisfy the relation (j43]) . the effective equation of state is 

V6. 



-1 



-Xxi 



(48) 



In what follows we shall discuss two fixed points separately. 
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1. Scaling solutions 

First we recall that Eq. (p6|) also gives xi = 2(q+a) • F^'o™ Eq. (|48| the effective equation of state is given by 



Wcff 



Q 



(49) 



which is independent of the form of p{(j),X). The accelerated expansion occurs for Wcs < —1/3, i.e., Q > A/2 or 
Q < -A. 

From Eq. P5|) we get 



X3 



2Q{Q + X)+3g{Y) 
8(4Q + A) 



where 



Y ^ — - 

Once the form of g{Y) is specified, one gets X3 and Y from Eqs. ((50|) and ([5T|) . This then determines 

3 . „ . „ 24 



■X3 . 



2(Q + A)' 



:(g + 2gi), a 



GB 



+ A 



X3 . 



(50) 



(51) 



(52) 



Let us consider an ordinary scalar field with an exponential potential, i.e., the model given in Eq. (|19p . In this case 
we get 



X3 



6Q(Q + A) + 9 - y^[6Q{Q + A) + 9]^ - 192ca(4Q + A)(Q + A)^ 



48(4g + A) 

We have chosen a minus sign in Eq. (ISS)) to recover 0:3 ^ as a ^ 0. When Q = this reduces 



X3 



16A 



and also we obtain 



2A2 



32caA3 



27(1 - y^l - 64caAV27) 



1 64 X3 
1 ~ -caA3 



If caA'^ is much smaller than unity, we find the following approximate fixed points: 



together with 



X2 



3 
A2 



3 2 2 

2a^' -3-gCaA , 



1 ° \3 

1 caA 

27 



o 16 X 



(53) 



(54) 



(55) 



(56) 



(57) 



which, in the limit a — > 0, recovers the scaling solution derived in Ref. [12l |3 

The above scaling solution is attractive to be used in a matter era to alleviate the coincide nt p roblem, but it does 
not give way to a late-time acceleration since the scaling solution is a global attractor [l^[l^[33|- If > A, however, 
the presence of the GB term can lead to the exit from the scaling matter era, as we will discuss later. 
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lig + 2gi) + 8^6x1X3 ^ 1 . (58) 



2. Scalar-field & GB dominated point 
Since fi^ + JIgb = 1 in this case, this gives 

X 

Making use of Eq. ((44|l together with this equation, we obtain 

16X3 = A[x?(5 + 2.gi) + 1] - 2V6{g + gi)xi . (59) 

If we specify the form of g(Y), one can get xi and x^ from the above equations. Let us consider the standard field 
with an exponential potential, i.e., p = 1 — c/Y . Then from Eq. (j59p we get 

A - VQxi 

X3 = ^ . (60) 

4(2 + V6Axi) 

Using Eq. ([55]) one obtains the following equation for xi: 

3(A - V6xi){V6Xxl - lOxj + \/6Aa;i - 2) + 4ca(2 + VGXxi)^ = . (61) 

When a = this has a solution xi — X/V6, which can lead to an accelerated expansion for A < ^/2. Let us obtain the 
solution for perturbatively under the assumption that a is much smaller than unity. Substituting xi = A/-\/6 + e 
for Eq. (|6ip . we find the following approximate relation 

A 4^/6ca(A2 + 2) /TT A^" 2ca 

^^ = 71" 3(6-A2) ' "^^V^l ^^J' "^ = 6^' 



and 



1 6cQ! A 

^^GB = g _ ^2 ' 1 - ^^GB ■ (63) 

Note that the effective equation of state is given by 

^ A2 8mA(A2 + 2) ^ A^ A2 + 2^ 

"- = -i + y-^(6^ = -^ + y-^"«- 

Thus the contribution of the GB term tends to reduce Wcs- In the absence of the GB coupling the late-time acceleration 
occurs for A^ < 2. In this case the scaling matter era is not present, since the existence of it requires the condition 
A^ > 3 [13 . Hence this case is not attractive to solve the coincident problem, although the standard matter-dominated 
epoch is followed by an accelerated expansion due to a shallow potential as in the case of a cosmological constant. 
When the GB term is present, there is a contribution of its energy density to the late-time acceleration. However, 
as we will see in a later section, the contribution of the GB term is restricted to be small from the constraint that it 
does not disturb the cosmological evolution during radiation and matter eras. 

B. Case of X fi: Exit from the scaling regime 

Strictly speaking, scaling solutions are present only for X — ^ m the presence of the GB coupling. Even when 
X ^ fi, however, they can exist approximately as long as the contribution of the GB term is negligibly small. This 
corresponds to a situation in which 2:3 is very much smaller than 1. Note that X3 can not be exactly zero since the 
relation between X2 and X3 is 

xlx3 = ^e'^^~'^^. (65) 

Still one can regard « as an approximate fixed point satisfying this relation. Then we approximately obtain 
the following 4 fixed points (A), (B), (C) and (D) using the results of Refs. [l^ [l^ recently obtained for X3 = 0. 
Hereafter, when we write X2 — or X3 — 0, it means that they are not exactly zero but are very small values. 
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1. K3 = and V&^xi/2 = -H/H'^ 



• (A) Scaling solutions. 
They are characterized by 



V6 Q Q{Q + A) + 3p,x 2 n 

^^=2(qTA)' ^°« = -0Ta' = (Q + A)2 ' 5(n = -3Q(Q + A), (66) 



where p_x = dp/dX. Note that xi and are the same as in the case where the GB term is present. This 
solution is stable when the following conditions are satisfied [l6j : 

^ <n^<l, A>0, (67) 



Q + \ 

where A is defined in Eq. (|3T|) . The latter condition automatically holds when we impose the stability of quantum 
fluctuations of the scalar field [lOl]. When Q — the scaling solution is stable for < $7^ < 1, as required for 
the existence of itself. 

• (B) Scalar-field dominated solutions. 
They are characterized by 

0,^ = 1, Wcff = W0 = -IH ■ (68) 

When Q is positive, this solution is stable for xi < \/6/[2{Q + A)]. 



2. X3 = and X2 = 

The solutions which satisfy X3 = and X2 = correspond to kinematic fixed points. These exist only when 
g(Y = XI/X2) is non-singular, namely, when g is expanded as 

9 = co-J2 cn^"" ' (69) 

n>0 

where cq and c„ are constants. Since gn{x2 ^ 0) = in this case, one can easily get the following points from 
Eq. jHSl). 

• (C) (/)MDE solution. 
This is characterized by 

n o 

Xi=-— , X2=0, ilci, = WcS = — . (70) 

6Cq 3co 

When Cq > 0, i.e., corresponding to a non-phantom field, this solution is a saddle point. Note that when Q — 
the (/)MDE is equivalent to a standard matter-dominated epoch. 

• (D) Pure kinetic solutions. 
These solutions are 

Xi^±l/y/c^, X2=0, fi0=Wcff = l- (71) 

We need positive cq for their existence. These are saddle or unstable nodes if Q > 0. Since fl^ — = 1, one 
can use pure kinetic solutions neither for matter/radiation eras nor for dark energy eras. 

From Eq. and ([1T|) we find that there are two other cases. 
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3. VGiixi/2 = -HjH^ and a;2 = 

• (E) Kinetic and GB dominated solutions. 

Since 2:2 = in this case, the form of the function g{Y) is restricted to be in the form ([69]) . When Q = and 
Co = these solutions satisfy 

F(xi) = 6M'xt-(24 + V6)Ma;? + [6(Ai2 + 5)4V6]a;2-5V6/iXi+6 = 0, X3 = ^^i(^-^^i) ^ (73) 



8(3^X1 - V6) ' 



together with an effective equation of state: 



V6 



-IJ.X1 



(73) 



We have w^s = — 1/3 for Xi = \/6/(3/i). Equation ((72|l possesses two real solutions. If 11 = 10, for example, 
we get {xi,X3, lUoff) = (0.08235, 1.40993, -0.3276), (0.12454, 0.02236, 0.01689). In both cases we do not have an 
accelerated expansion {wcS < —1/3). When ^ > we find that the values of xi corresponding to two real 
solutions of Eq. ((72|l are larger than V6/(3/i), which means Wcg > — 1/3 from Eq. ([73|) . Hence one can not use 
these solutions for dark energy. 

These fixed points correspond to the absence of the field potential, in which case the accelerated solution has 
not been found in Ref. [23 • This is consistent with the above result. 



4. = -H/H^ and V6\xi/2 = -H/H^ 



(F) de-Sitter point. 



When X fi there exists a fixed point xi = from Eqs. (f40|) and (f4T|) . Since H/H^ = in this case, this corresponds 
to a de-Sitter solution. In fact one has WcS = — 1 from Eq. (I48p . Equation (|38p gives 



(3 + 3gxl) + I2AX3 ~V6{g + gi)xi + {Q + X)ig + 2gi)xl - Sx^, - 



= 0. 



Meanwhile Eq. (139]) leads to 



A 



-V6(g + gi)xi + (Q + X){g + 2gi)xl - 8x3 - 



= 0. 



Then we find that the fixed point satisfies 

gxl = -1 , .Ti = 0. 

Once the form of g{Y) is specified, one can get X2 by solving this equation. From Eq. (|75p we get 

1 r. 



2:3 



{Q + X){g + 2gi)xi 



+ 9i)xi 



(74) 



(75) 



(76) 



(77) 



Let us consider a standard scalar field with an exponential potential, i.e., the model given by Then the above 
de-Sitter solution corresponds to 

(a;i,a;2,X3) = (0,l/VS,A/8). (78) 
More generally Eq. ([76]) can be satisfied when g = g{Y) is written in the form 

g{Y)^[co-{c/Y)]f{Y), (79) 
where f{Y) is a function which approaches a constant d as xi — > 0. In this case we obtain the following fixed point: 

{xuX2,X3) = {0,l/V^,X/8). (80) 
The model (|79p includes a tachyon field with a Lagrangian given in Eq. (I2ip . 



11 



In order for the existence of the de-Sitter solution, it is crucial to have the term cjY in the expression of For 
example, when g{y^ is written as a sum of positive powers of such as a dilatonic ghost condensate model given in 
Eq. ((20|) . we do not have such a de-Sitter solution. 

From Eqs. ([55]) and ([50)1 we find that this fixed point satisfies 

1]^ = 1, 1^GB=0. (81) 

Since the field is frozen at the de-Sitter point, the GB term does not have any contribution to the energy density of 
the universe. 

The appearance of the de-Sitter point comes from the presence of the GB term. For example when g{Y) — 1 — c/Y , 
Eq. ([50]) reduces to 

+ iH^ - cXer^'^ + 2A f{(j}){H'^ + ij) = , (82) 

where we dropped the couphng Q. Since f'{(p)H^ — X/8 and iJc = at the critical point, the effective potential 
Vcs{4>) for the field satisfies the relation 

^ = -cAe-^^ + 3Ail2 . (83) 
dcj) 

The last term appears because of the presence of the GB term, which gives rise to a potential minimum for the field 
(j). Then after the field drops down at this minimum, the universe exhibits a de-Sitter expansion. 

Let us consider the stability of the de-Sitter solution. By perturbating Eqs. (IHS)) . and (|¥T|) about the fixed 
point, we obtain a 3 x 3 Jacobian matrix M for perturbations Sxi, 6x2 and 6x3 (we neglect the contribution of 
radiation). For the model (flQl) the eigenvalues of the matrix Ai are 



Ai — —3, A2,3 — 2 



-1 ± 



8X{X- n) 



(84) 



3(2 + 3A2) 

This explicitly shows the following property for the stability of the de-Sitter point: 

• (i) Stable for /i > A. 

• (ii) Saddle for /i < A. 

When fj, = X one has Ai = —3, A2 = and A3 = —3, which means that the de-Sitter solution is marginally stable. 

The same eigenvalues as given in Eq. (|84p can be obtained by demanding the conditions fl^ = 1, {g -\- gi)xi — 
and A — 1 a.t the fixed point without specifying the form oi g{Y). This includes the tachyon model pT|) with a choice 
c = 1. Provided that /i > A the system falls down to the stable de-Sitter point. 



IV. COSMOLOGICAL DYNAMICS 



In this section we shall study cosmological dynamics for the scaling Lagrangian (fT5|) with the GB coupling given 
in Eq. psp . Our interest is to find a case in which a late-time accelerated expansion is realized by the presence of 
the GB term. Such a possibility can be accommodated by using the de-Sitter fixed point discussed in the previous 
section. We would also like to make use of scaling solutions during the matter-dominated era in order to alleviate the 
coincident problem. 

When A = /i there exists a scaling matter era with < $7^, ficB < 1, but the system does not get away from the 
scaling regime to give rise to a late-time acceleration. We require that fi does not equal to A in order to exit from the 
scaling matter era. Let us consider the approximate scaling solution (j66p which exists under the conditions 0:3 ~ 
and •\/6Axi/2 = —H/H^. Perturbing Eq. (j4T|l about the fixed point, we obtain 

To get Wcs — during the matter era, we require Q <^ X from Eq. (p^ . Since we are now considering positive A, 
the system departs from 0:3 = for ^ > A whereas it approaches x^ = Q for ji < X. As we mentioned in the previous 
section, the stability of the perturbations 6x1 and 6x2 is ensured when the condition —Q/{Q -I- A) < fi^ < 1 is 
satisfied, i.e., A(A -I- Q) > Sp x > —2Q{X + Q). For a non-phantom field {p x > 0) with Q = 0, this condition reduces 
to A^ > 3p.x- From the above discussion we find the following property for the stability of the scaling solution: 
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FIG. 1: Evolution of the variables Q,^, J^gb, fJm and Q,r together with the effective equation of state Wcff for A = 4, ^ = 12, 
a — 10"'^^, (5 = and c = 10~^ in the case of an ordinary field with an exponential potential. We choose initial conditions 
xi — 10~*, X2 = 10~^, = 3.6 X 10"'' and X4 = 0.999. The solution is in a scaling regime during the matter-dominated epoch 
and finally approaches the de-Sitter universe characterized by Q.^ = 1, Qgb ~ and uieff = — 1. The energy fraction of the field 
(/!> during the scaling regime is = S/A^ = 0.1875. 



• (i) Saddle for ^ > A and A(A + Q) > 3p.x > -2Q(A + Q). 

• (ii) Stable for ^ < A and A(A + Q) > 3p^x > -2Q{\ + Q). 

Then in the case (i) the solutions can exit from the scaling regime to connect to the dark energy era. 

As we showed in the previous section, the de-Sitter point (F) is stable for > A, whereas it is a saddle point for 
jjL < \. Then it is clear that the scaling solution can be connected to the de-Sitter solution provided n > X and 
A(A + Q) > 'ip,x > —2Q{X + Q). The de-Sitter solution is present only for a restricted class of models. In what 
follows we shall investigate cosmological evolution for such two models: (i) an ordinary field and (ii) a tachyon field. 



A. An ordinary field with an exponential potential 

We shall first study the case of an ordinary field with an exponential potential {p = X — ce~^'^). Since p,x = 1 in 
this case, one can have a scaling solution in the matter era if A^ > 3 in the absence of the coupling Q. In Fig. [1] we 
plot the evolution of fi^, flcB, and fi^ together with WeS for X — 4, — 12, a = 10~^^ and Q — 0. The system 
starts from a radiation-dominated epoch, which is followed by the scaling regime in the matter-dominated era. Since 
/i > A the solution exits from the scaling regime and finally approaches the stable de-Sitter fixed point. The final 
attractor solution actually satisfies fl^ = 1 and ^gb = 0, as estimated analytically. 

Figure [T] shows that the energy fraction of the GB term grows right after end of the matter era, which begins to 
decrease after the increase of ft^. In Fig. [T]the effective equation of state temporally takes a local minimum value 
WcS ~ —0.6 around Qm ~ 0.3-0.4, which means that the accelerated expansion indeed occurs at the present epoch in 
this scenario. Note that the rapid transition of Wcs just after the matter-dominated era is associated with the growth 
of r^GB- As /i is increased, WcS tends to get smaller. When A = 4, for example, we find that the phantom equation of 
state, WcS < —1, is realized for /i > 25. Meanwhile the increase of /i leads to a shorter period of the matter-dominated 
epoch. Hence the transient phantom stage is realized at the expense of such a short matter period. If we take larger 
A, it is also difficult to get smaller values of WeS satisfying the condition for an accelerated expansion compatible with 
observations. Thus A is bounded from above as well. 

In Refs. [1^ Koivisto and Mota placed observational constraints using the Gold data set of Supernova la [3| 
together with the CMB shift parameter data [i^. The parameter A is constrained to be 3.5 A < 4.5 at the 90% 
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FIG. 2: Evolution of the same variables as in Fig.[T]but for A = 1, ^ = 10^, a = 10"^^, Q = and c = 10"^. We choose initial 
conditions xi = 10~^^, X2 = 10"^'', x-j, — 3.3 x 10"'* and xt, = 0.999. After a standard matter era without a scaling regime, 
the system falls down to a de-Sitter fixed point before it reaches a scalar-field dominated point. The efi'ective equation of state 
temporally drops down toward the phantom region {waS < — 1). 

confidence level, see Fig. 3 in Ref. ^2^^. If the solutions are in the scaling regime in radiation era the constraint 
coming from Nucleosynthesis gives A > 4.47 [1, [i3|i in which case the allowed range of A is severely constrained. 

When A = we have shown that there exists a scalar- field & GB dominated point with fi^ -t- r^GB = 1- This 
solution can be also used for a late-time acceleration. Note that in the absence of the GB term the existence of the 
stable scalar-field dominated solution demands the condition A^ < 3 [l^l, in which case the scaling solution with 
< 1 does not exist in the matter-dominated epoch. To get a stable scalar-field dominated solution, one has to use 
the (/)MDE solution ([70]) in the preceding matter period (the standard matter era corresponds to Q = 0). In this case 
X2 and X3 need to be very small during a radiation era to have Jl^ <C 1 and SIgb ^ 1, which restricts the coupling 
a very small through the relation x^x^, = a/3. Then, at the scalar-field & GB dominated point, Ogb in Eq. is 
negligibly small relative to fi^, in which case the effect of the GB term is not important around the present epoch. 
Thus, for small A, the 0MDE is followed by the scalar-field & GB dominated point, but the impact of the GB term 
for the late-time acceleration is not significant for A = 

Then what happens for fi > X 7 In this case the GB term finally comes out to lead to the de-Sitter expansion 
by giving rise to a minimum of an effective potential. If fi is not much different from A, this occurs at sufficiently 
late-times after the solutions approach the scalar-field dominated point. On the other hand, if /i » A, the system 
can be trapped at the de-Sitter point before it reaches the scalar-field dominated solution. In Fig. [2] we plot the 
cosmological evolution for X — 1, fi — 100, a — 10^^^ and Q = 0. In this case the scaling solution does not exist in 
the matter era because of the smallness of A (< \/3). We find that the GB term becomes important around fl^ = 0.7, 
which leads to the rapid decrease of WeS toward the phantom region. The effective equation of state approaches — 1 
after the field settles down at the potential minimum. The effect of the GB term can be seen at present time in such 
a case, but the absence of the scaling matter era is not attractive to alleviate the coincidence problem. 

B. A tachyon field with an inverse square potential 

As we have already mentioned, the tachyon field also possesses the de-Sitter fixed point. Meanwhile, when Q — 0, 
the scaling solution does not exist if the background fluid corresponds to a non-relativistic matter {wm = 0). In fact. 




^ Note that the definition of A in Ref. [26| is different from ours by the factor \/2. 
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for the function g{Y) given in Eq. (|2ip. Eq. (l66l) shows that the scahng solution satisfies 

3c 



(Q + A)2 



QiQ + \) 



(86) 



^\/r^2F=^0(Q + A). (87) 

When (3 = one gets Y ^ 1/2 from Eq. (|87p. which then means that il^ diverges from Eg. ( [SS]). More generaUy we 
have flrj, = 3{1 + Wm) / \/ —Wm) for the background fluid with an equation of state Wm [4?, '4^, showing that the 
scahng solution exists only for Wm < 0. Hence in the absence of the coupling Q one can not have a successful scaling 
epoch followed by a de-Sitter expansion. 

On the other hand the scaling matter era can exist if the coupling Q is present. As we already mentioned, we require 
Q <C A to recover the equation of state WcS ~ 0- In this case we get the approximate relation ^1 — 2Y « QX/3c from 
Eq. ([87|) . Then from Eq. ([86]) wc find that fi^ is approximately given by 

o 9c2 + (gA)2 



QA3 

For example, when Q = 0.1, A = 10 and c = 1, one has = 0.1 and WcS = —0.0099. Note that the variable A 
defined in Eq. (pij) is given hy A = (1 — 2y)^/^/c, which is positive for c > (corresponding to a positive potential). 
Then the scaling solution is stable provided that —Q/{Q + A) < il^ < 1 and c > 0. It becomes a saddle point if the 
GB coupling is present with /i larger than A. In this case the scaling matter era is followed by the de-Sitter solution. 

In Fig. [3] we plot the cosmological evolution for Q = 0.1, X — 7, fj, — 20 and a = lO^'^^. We find that the solution 
first enters the scaling matter era during which Cl„i and fl^ are constants. This regime is followed by the growth of 
the GB term, which leads to a rapid decrease of Wcs toward WcS ~ —0.7 around the present time (f2„i — 0.3-0.4). 
Numerically we checked that the minimum value of WcS temporarily reached after the scaling matter era gets smaller 
if we choose smaller A or larger fi. Finally the solution falls down to the de-Sitter fixed point with — 1, CIqb — 
and WcB = — 1- 

From the above discussion we require that Q is bounded from below to get the scaling matter era and that it is 
bounded from above to obtain the effective equation of state which does not differ from WeS = much. It is certainly 
of interest to find the parameter spaces of Q and A to realize the sequence of the scaling matter era and the de-Sitter 
expansion satisfying observational constraints, which we leave for future work. 



Ghost conditions 



Recently a number of authors [3l|, |32| discussed ghost conditions in GB cosmologies. It is possible to see the 
signature of ghosts by studying gravitational perturbations. We shall consider scalar and tensor perturbations about 
the FRW background [i^: 



-(1 + 2A)dt'^ + 2ad,Bdx'dt + c? [(1 + 2?A)% + 2d^jE + 2h,j] dx'dx^ , 



(89) 



where di represents the spatial partial derivative d/dx^. Here A, B, and E are scalar quantities, whereas hij is a 
tensor quantity. We shall investigate whether or not ghosts appear for an ordinary scalar field with an exponential 
potential, i.e., p — X — ce"^'^. In what follows we neglect the contribution of the matter fluid. Although the equation 
of scalar perturbations is modified in the presence of matter, the equation of tensor perturbations is unchanged. 

Defi ning the gauge-invariant comoving perturbation, TZ = ^ — H6(f>/(j), the Fourier modes of scalar perturbations 
satisfy HI mi 



1 P 

(a^QsTZy +4—71 = 0, 



where A: is a comoving wavenumber and 

(1 - 8Hf) \{1 - 8i/ /)(/)2 + m{H^ff 



i/2(i - l2HfY 



2 



(1 - 8Hf) 


(1 - 8i//)02 + 96(i?V)^ + 12SH^Hf^ 


+ 256iH^f)\f~Hf) 


(1 - 8Hf) 


(1 - 8Hf )<j)^ + 96(i?2/)2 





(90) 

(91) 
(92) 
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FIG. 3: Evolution of the variables fl^, Qgb, fim and Qr together with the effective equation of state uicff for Q = 0.1, A = 7, 
H — 20, a — 10~^^ and c = 1 in the case of a tachyon field with an inverse square potential. Initial conditions are chosen to 
be xi = 10~®, X2 = 10"'*, X3 = 3.0 x 10~^^ and X4 = 0.999. The scaling matter era is present because of the presence of the 
coupling Q. 



Decomposing tensor perturbations into eigenmodes of the spatial Lagrangian, V^e^ = —k'^eij, with scalar amplitude 
h{t), i.e., hij = h(t)eii, where have two polarization states, the Fourier modes of tensor perturbations obey the 
equation of motion [43, 



^2 

{a^QThy +4^h = 0, (93) 



where 



Qt = 1-8H/, c| = -1-^. (94) 

The perturbations exhibit exponential instabilities if Cg and are negative. The propagation speeds of scalar and 
tensor modes are superhiminal if Cg and are larger than unity. In Refs. [sil Is^ it was shown that the no-ghost 
state to ensure a consistent quantum field theory gives tlic constraints 1 — 8Hf > and 1 — 8/ > 0. Thus stable, 
non-superluminal and no ghost states require tlic conditions 

0<c^<l, 0<c|<l, e = l- 8i// > . (95) 

In Fig. [5] we plot the evolution of c| and for the model p = X ~ ce^^'^ with the same model parameters and 
initial conditions as given in Fig. [TJ During the scaling matter era in which is much smaller than of order unity 
[see Eq. ([56]) ]. we have that \Hf\ <C 1 and |/| <^ 1 giving Cg ~ 1 and ~ 1. For the de-Sitter point ([78]) one can 
easily show that Cg = = 1. These properties are in fact confirmed in the numerical simulation of Fig. [J] 

We find from Fig. [¥] that Cg and temporally become negative during the transition from the scaling matter 
era to the final de-Sitter era. This corresponds to the stage in which the contribution of the GB term is dominant. 
The perturbations exhibit exponential instabilities associated with the appearance of ghosts. The propagation speed 
of scalar perturbations can be modified in the presence of a matter fluid, but that of tensor perturbations is not 
affected. Hence the negative instability of tensor perturbations is robust even taking into account the contribution of 
the matter. This is consistent with the results of Ref. ^] in which the negative instability of tensor perturbations 
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FIG. 4: Evolution of the variables ^ = 1 — 8Hf, Cg and Cx for an ordinary field with an exponential potential with the same 
model parameters and initial conditions as given in Fig. [T] This shows that both Cg and are temporally negative around the 
transition from the scaling matter era to the final de-Sitter era. Moreover the propagation speed of the scalar mode becomes 
superluminal. 

has been also found if inflationary solutions are constructed by the dominance of the GB term. In Fig. 3] we also find 
that the propagation speed of scalar perturbations temporally become superluminal^ . 

For the model parameters given in Figs. [1] and [4] the propagation speed ct becomes imaginary before the system 
reaches the present epoch (while Cg is positive). Meanwhile if we choose smaller values of /i, it is possible to avoid 
reaching the negative values of c^. For smaller n, however, the minimum value of WcS tends to be larger. Hence the 
exponential instability of tensor perturbations can be avoided at the expense of losing rapid accelerated expansion 
of the universe. It is of interest to find parameter ranges of A and /i which lead to positive and also satisfy the 
observational constraints of WcS- 

V. CONCLUSIONS 

In this paper we have discussed the viability of Gauss-Bonnet (GB) dark energy models which possess cosmological 
scaling solutions. In order to alleviate the coincident problem of dark energy, it is of interest to find a cosmological 
trajectory which is in a scaling regime (fl^/flm ^constant) during a matter era and finally exits to give rise to a 
late-time accelerated expansion. We have tried to understand the property of scalar-field dark energy models which 
allow such a possibility. 

First we derived the form of the GB coupling together with the form of the scalar-field Lagrangian for the existence 
of scaling solutions by starting from a very general action given in Eq. ([1]). The form of the GB coupling has been 
found to be /(0) c>c e^'^ together with the field Lagrangian density p — Xg{Xe^'^), where A is defined in Eq. (fT6|) and 
g is an arbitrary function. The field Lagrangian density takes the same form as in the case where the GB coupling is 
absent. 

We have also derived autonomous equations and fixed points for the scaling Lagrangian p — XglXe^"^) with the 
GB coupling /((/>) c>c e^"^. When A = /i, we obtained the scaling solution characterised by xi = \/&/[2{Q + A)] and 




^ In the context of K-essence, Bonvin et al. [49(| recently found that the propagation speed of scalar perturbations also temporally becomes 
superluminal in order to solve the coincident problem of dark energy. 
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the scalar-field dominated solution characterised by fl^f, + floB = 1 without specifying any form of the Lagrangian. 
In this case, however, the solutions do not get away from the scaling matter regime to connect to a dark energy era. 
When /i > A the presence of the GB term can lead to the exit from the scaling regime. 

The GB term allows an interesting possibility to have a late-time de-Sitter fixed point. This point is present if the 
solutions for gx^ — and xi = exist. The models which possess this de-Sitter solution are characterised by the form 
g{Y) = {cq — c/Y)f{Y), where f{Y) is a function that approaches a constant as xi 0. The ordinary field with an 
exponential potential and the tachyon field with an inverse square potential (pij) belong to this class, whereas 
the dilatonic ghost condensate model (PO)) does not. 

We have studied cosmological evolution for the above two models in which the de-Sitter point is present. In the 
case of an ordinary field with an exponential potential one can in fact obtain a viable trajectory along which the 
scaling matter period is followed by the de-Sitter solution provided that /i > A and A^ > 3 for Q = 0. There is also 
an interesting possibility in which the effective equation of state WcS temporarily becomes smaller than —1 because of 
the growth of the GB term. This is realized by choosing larger values of /i, but it also corresponds to a shorter matter- 
dominated period. The late-time acceleration is possible even for A = ^ by using a scalar-field & GB dominated point 
{ilcj, + floB = 1), but in this case the contribution of the GB term is very small at the present epoch. Moreover we do 
not have a scaling matter era in this case, which is not a welcome feature in solving the coincidence problem. 

The tachyon field has a de-Sitter fixed point, but in this case the scaling solution is absent for a non-relativistic 
background fluid {wm = 0) if Q = 0. Hence we do not have a successful sequence of the scaling matter era and the 
de-Sitter expansion. If the coupling Q is present, however, it is possible to realize this sequence when fi is larger than 
A. The values of Q are bounded from both above and below to get a scaling matter era with WcS close to 0. An 
alternative approach to the problem discussed here can be provided by Renormalization Group (RG) equations. The 
scaling solution presented here turns out to be the fixed point of RG flow equations; the form of the field Lagrangian 
and the GB coupling are automatically fixed in this method [soj . 

We also discussed ghost conditions for the GB dark energy models by considering gravitational perturbations. For 
the ordinary field with an exponential potential we found that the propagation speeds of scalar and tensor modes 
can be imaginary during the transition from the scaling matter era to the final de-Sitter attractor. This is associated 
with the appearance of ghosts which lead to negative instabilities of perturbations. This ghost stage can be avoided 
by tuning the model parameters at the expense of obtaining a minimum effective equation of state closer to —1. 

In this paper we have not studied local gravity constraints on the GB coupling fsT', 'E^ . If the contribution of the GB 
coupling is large at present epoch, this may contradict with local gravity experiments unless an accidental cancellation 
occurs to avoid the large variation of effective gravitational constant. We leave future work for the construction of 
scaling GB dark energy models consistent with such constraints. We hope that the results obtained in this paper will 
be useful to find such viable models. 



ACKNOWLEDGEMENTS 



We thank L. Amendola, S. Jhingan, S. Nojiri and S. D. Odintsov for useful discussions. M. S. is supported by ICTP 
and lUCAA through their associateship programs. S.T. is supported by JSPS (Grant No. 30318802). 



[1] S. Perlmutter et al, Astrophys. J. 517, 565 (1999); A. G. Riess et ai, Astron. J. 116, 1009 (1998); Astron. J. 117, 707 
(1999); J. L. Tonry et al, Astrophys. J. 594, 1 (2003); R. A. Knop et al, Astrophys. J. 598, 102 (2003). 

[2] V. Sahni and A. A. Starobinsky, Int. J. Mod. Phys. D 9, 373 (2000); V. Sahni, Lect. Notes Phys. 653, 141 (2004); 
S. M. Carroll, Living Rev. Rel. 4, 1 (2001); T. Padmanabhan, Phys. Rept. 380, 235 (2003); P. J. E. Peebles and 
B. Ratra, Rev. Mod. Phys. 75, 559 (2003). 

[3] E. J. Copeland, M. Sami and S. Tsujikawa, arXiv:hep-th/0603057. 

[4] S. Weinberg, Rev. Mod. Phys. 61, 1 (1989). 

[5] Y. Fujii, Phys. Rev. D 26, 2580 (1982); L. H. Ford, Phys. Rev. D 35, 2339 (1987); C. Wetterich, Nucl. Phys B. 302, 668 
(1988); B. Ratra and J. Peebles, Phys. Rev D 37, 321 (1988); Y. Fujii and T. Nishioka, Phys. Rev. D 42, 361 (1990); 
T. Chiba, N. Sugiyama and T. Nakamura, Mon. Not. Roy. Astron. Soc. 289, L5 (1997); R. R. Caldwell, R. Dave and 
P. J. Steinhardt, Phys. Rev. Lett. 80, 1582 (1998); I. Zlatev, L. M. Wang and P. J. Steinhardt, Phys. Rev. Lett. 82, 896 
(1999); B. Boisseau, G. Esposito-Farese, D. Polarski and A. A. Starobinsky, Phys. Rev. Lett. 85, 2236 (2000). 

[6] T. Chiba, T. Okabe and M. Yamaguchi, Phys. Rev. D 62, 023511 (2000); C. Armendariz-Picon, V. Mukhanov, and P. J. 
Steinhardt, Phys. Rev. Lett. 85, 4438 (2000); Phys. Rev. D 63, 103510 (2001). 

[7] A. Sen, JHEP 9910, 008 (1999); JHEP 0204, 048 (2002); M. R. Garousi, Nucl. Phys. B584, 284 (2000); Nucl. Phys. 
B 647, 117 (2002); G. W. Gibbons, Phys. Lett. B 537, 1 (2002); T. Padmanabhan, Phys. Rev. D 66, 021301 (2002); 



18 



T. Padmanabhan and T. R. Choudhury, Phys. Rev. D 66, 081301 (2002); M. R. Garousi, M. Sami and S. Tsujikawa, Phys. 

Rev. D 70, 043536 (2004). 
[8] M. Gaspcrini, F. Piazza and G. Vcneziano, Phys. Rev. D 65, 023508 (2002). 
[9] N. Arkani-Hamed, P. Creminelli, S. Mukohyama and M. Zaldarriaga, JCAP 0404, 001 (2004). 
[10] F. Piazza and S. Tsujikawa, JCAP 0407, 004 (2004). 

[11] P. G. Fcrrcira and M. Joyce, Phys. Rev. Lett. 79, 4740 (1997); Phys. Rev. D 58, 023503 (1998). 
[12] E. J. Copcland, A. R. Liddle and D. Wands, Phys. Rev. D 57, 4686 (1998). 

[13] A. P. Billyard, A. A. Coley and R. J. van den Hoogen, Phys. Rev. D 58, 123501 (1998); A. R. Liddle and R. J. Scherrer, 
Phys. Rev. D 59, 023509 (1999); J. P. Uzan, Phys. Rev. D 59, 123510 (1999); R. J. van den Hoogen, A. A. Coley and 
D. Wands, Class. Quant. Grav. 16, 1843 (1999); A. Nuncs and J. P. Mimoso, Phys. Lett. B 488, 423 (2000); A. dc la 
Macorra and G. Piccinelh, Phys. Rev. D 61, 123503 (2000); S. C. C. Ng, N. J. Nunes and F. Rosati, Phys. Rev. D 64, 
083510 (2001); C. Rubano and J. D. Barrow, Phys. Rev. D 64, 127301 (2001); Z. K. Guo, Y. S. Piao, R. G. Cai and 
Y. Z. Zhang, Phys. Lett. B 576, 12 (2003); Z. K. Guo, Y. S. Piao and Y. Z. Zhang, Phys. Lett. B 568, 1 (2003); S. Mizuno, 
S. J. Lee and E. J. Copeland, Phys. Rev. D 70, 043525 (2004); E. J. Copeland, S. J. Lee, J. E. Lidsey and S. Mizuno, Phys. 
Rev. D 71, 023526 (2005); M. Sami, N. Savchenko and A. Toporensky Phys. Rev. D 70, 123526 (2004); V. Pettorino, 
C. Baccigalupi and F. Perrotta, JCAP 0512, 003 (2005); L. Amendola, S. Tsujikawa and M. Sami, Phys. Lett. B 632, 155 
(2006); Y. Gong, A. Wang and Y. Z. Zhang, Phys. Lett. B 636, 286 (2006). 

[14] S. Tsujikawa and M. Sami, Phys. Lett. B 603, 113 (2004). 

[15] S. Tsujikawa, Phys. Rev. D 73, 103504 (2006). 

[16] L. Amendola, M. Quartin, S. Tsujikawa and I. Waga, Phys. Rev. D 74, 023525 (2006). 
[17] T. Barreiro, E. J. Copeland and N. J. Nunes, Phys. Rev. D 61, 127301 (2000). 
[18] V. Sahni and L. M. Wang, Phys. Rev. D 62, 103517 (2000). 
[19] A. Albrecht and C. Skordis, Phys. Rev. Lett. 84, 2076 (2000). 

[20] A. A. Coley and R. J. van den Hoogen, Phys. Rev. D 62, 023517 (2000); S. A. Kim, A. R. Liddle and S. Tsujikawa, Phys. 

Rev. D 72, 043506 (2005). 
[21] A. R. Liddle, A. Mazumdar and F. E. Schunck, Phys. Rev. D 58, 061301 (1998). 

[22] J. E. Lidsey D. Wands and E. J. Copeland, Phys. Rept. 337, 343 (2000); M. Gasperini and G. Veneziano, Phys. Rept. 
373, 1 (2003) 

[23] I. Antoniadis, J. Rizos and K. Tamvakis, Nucl. Phys. B 415, 497 (1994); M. Gasperini, M. Maggiore and G. Veneziano, 
Nucl. Phys. B 494, 315 (1997); S. Foffa, M. Maggiore and R. Sturani, Nucl. Phys. B 552, 395 (1999); N. E. Mavromatos 
and J. Rizos, Phys. Rev. D 62, 124004 (2000); JHEP 0207, 045 (2002); C. Cartier, E. J. Copeland and R. Madden, JHEP 
0001, 035 (2000); S. Tsujikawa, R. Brandenberger and F. Finelli, Phys. Rev. D 66, 083513 (2002); S. Tsujikawa, Class. 
Quant. Grav. 20, 1991 (2003). 

[24] E. Elizalde, S. Nojiri and S. D. Odintsov, Phys. Rev. D 70, 043539 (2004); S. Nojiri and S. D. Odiutsov, Phys. Lett. B 
595, 1 (2004); Phys. Rev. D 70, 103522 (2004); S. Nojiri, S. D. Odintsov and S. Tsujikawa, Phys. Rev. D 71, 063004 
(2005). 

[25] S. Nojiri, S. D. Odintsov and M. Sasaki, Phys. Rev. D 71, 123509 (2005). 
[26] T. Koivisto and D. F. Mota, arXiv:astro-ph/0606078. 

[27] G. Calcagni, S. Tsujikawa and M. Sami, Class. Quant. Grav. 22, 3977 (2005); M. Sami, A. Toporensky, P. V. Tretjakov 

and S. Tsujikawa, Phys. Lett. B 619, 193 (2005). 
[28] B. M. N. Carter and I. P. Neupane, Phys. Lett. B 638, 94 (2006); JCAP 0606, 004 (2006); I. P. Neupane, arXivihep- 

th/0602097. 

[29] S. Nojiri and S. D. Odintsov, Phys. Lett. B 631, 1 (2005); arXiv:hep-th/0601213. 

[30] G. Cognola, E. Ehzalde, S. Nojiri, S. D. Odintsov and S. Zerbini, Phys. Rev. D 73, 084007 (2006). 

[31] A. De Felice, M. Hindmarsh and M. Trodden, JCAP 0608, 005 (2006). 

[32] G. Calcagni, B. de Carlos and A. De Felice, Nucl. Phys. B 752, 404 (2006). 

[33] S. Nojiri, S. D. Odintsov and M. Sami, Phys. Rev. D 74, 046004 (2006). 

[34] S. Tsujikawa, Annalen Phys. 15, 302 (2006) [arXiv:hep-th/0606040]. 

[35] A. K. Sanyal, arXiv:astro-ph/0608104. 

[36] K. Andrew, B. Bolen and C. A. Middleton, arXiv:hep-th/0608127. 

[37] L. Amendola, Phys. Rev. D 62, 043511 (2000). 
[38] L. Amendola, Phys. Rev. D 60, 043501 (1999). 

[39] B. Gumjudpai, T. Naskar, M. Sami and S. Tsujikawa, JCAP 0506, 007 (2005). 

[40] L. R. W. Abramo and F. FineUi, Phys. Lett. B 575, 165 (2003). 

[41] J. M. Aguirregabiria and R. Lazkoz, Phys. Rev. D 69, 123502 (2004). 

[42] E. J. Copeland, M. R. Garousi, M. Sami and S. Tsujikawa, Phys. Rev. D 71, 043003 (2005). 
[43] T. Koivisto and D. F. Mota, arXiv:hep-th/0609155. 

[44] A. G. Riess et al. [Supernova Search Team Collaboration], Astrophys. J. 607, 665 (2004). 

[45] C. J. Odman, A. Melchiorri, M. P. Hobson and A. N. Lasenby Phys. Rev. D 67, 083511 (2003). 

[46] H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 78, 1 (1984); V. F. Mukhanov, H. A. Feldman and R. H. Branden- 
berger, Phys. Rept. 215, 203 (1992); B. A. Bassett, S. Tsujikawa and D. Wands, Rev. Mod. Phys. 78, 537 (2006). 
[47] J. c. Hwang and H. Noh, Phys. Rev. D 71, 063536 (2005). 
[48] Z. K. Guo, N. Ohta and S. Tsujikawa, arXiv:hep-th/0610336. 
[49] C. Bonvin, C. Caprini and R. Durrer, Phys. Rev. Lett. 97, 081303 (2006). 



19 



[50] S. Jhingan et. al, work in progress. 

[51] L. Amendola, C. Charmousis and S. C. Davis, arXiv:hep-th/0506137. 
[52] L. Amendola, C. Charmousis and S. C. Davis, work in progress. 



